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Abstract. We introduce a family of real random polynomials of degree n 
whose coefficients are symmetric independent Gaussian variables with variance 
(a 2 ,) = e~ k , indexed by a real a > 0. We compute exactly the mean number 
of real roots (N n ) for large n. As a is varied, one finds three different phases. 
First, for < a < 1, one finds that (N n ) — (^logn. For 1 < a < 2, there is 
an intermediate phase where (N n ) grows algebraically with a continuously varying 



exponent, (N n ) ~ f y 9L ^- n a l 2 . And finally for a > 2, one finds a third phase where 
(N n ) ~ n. This family of real random polynomials thus exhibits a condensation of 
their roots on the real line in the sense that, for large n, a finite fraction of their roots 
(N n ) /n are real. This condensation occurs via a localization of the real roots around 
the values ±exp [f (k + \) a ~ l ], K k < n. 
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1. Introduction 

Since the early work of Bloch and Polya pQ in the 30's, the study of random algebraic 
equations has now a long story [2, [3]. In the last few years, it attracted a renewed 
interest in the context of probability and number theory |4j, as well as in the field of 
quantum chaos [5j. Recently, we showed that there are also interesting connections 
between random polynomials and persistence properties of physical systems [6], [7]. 

Here we consider real random polynomials, i.e. polynomials with real random 
coefficients, of degree n. While these polynomials have exactly n roots in the complex 
plane, the number of roots on the real line N n is a random variable. One would like to 
characterize the statistics of this random variable and a natural question is thus : what is 
the mean number (N n ) of real roots and how does it behave with n for large n [4]? This 
question has been widely studied in the past for Kac's polynomials K n (x) = J2k=o a k xk 
where a& are independent and identically distributed (i.i.d.) random variables of finite 
variance (a 2 .) = a 2 . In that case it is well known that (N n ) ~ - logn, independently of a. 
This result was first obtained by Kac [8] for Gaussian random variables and it was later 
shown to hold also for a wider class of distributions of the coefficients [21 E]. Interesting 
generalizations of Kac's polynomials have been studied in the literature where are 
independent Gaussian variables but non identical, such that (a|) = fc d ~ 2 , where d > is 
a real number, leading to (N n ) ~ 7r _1 (l + yd]2) logn [HE]- Given the robustness of this 
asymptotic logarithmic behavior of (N n ), it is natural to search for random polynomials 
for which (N n ) increases faster than logn, for instance algebraically. 

One such instance is provided by the real Weyl polynomials W n (x) defined by 

n u 

^n(*)=]Te ^ , (1) 

where €k are i.i.d. random variables of zero mean and unit variance. Thus here, 
a-k = Cfe/VH and the variance is (a 2 ,) = l/k\, which for large k behaves as (a 2 ,) oc e~ hlogk . 
For these real polynomials in Eq. (pQ), it is known that (N n ) oc n 1//2 . For instance, 
in the special case where e& are Gaussian random variables of unit variance, one has 
(N n ) ~ - y/n [Zl [JO]. Another interesting and intriguing instance of real random 
polynomials was introduced a long time ago by Littlewood and Offord [TT] who studied 
the random polynomials L n {x) given by 

L ^) = ^ + E e *(|rp ( 2 ) 

where = ±1 with equal probability. Thus in this case = ek/{k\) k and the variance 
is (a|) = l/(k\) 2k , which behaves for large k as (a\) oc e - 2fc2l °s fc . Using algebraic 
methods, they showed that such polynomials L n (x) have all their roots real and therefore 
(N n ) = n. 

We thus have here two examples of real random polynomials in Eq. (0Q) and Eq. ((2j) 
where, at variance with Kac's polynomials, (N n ) grows algebraically with n. In the 
second example (j2J), the number of real roots is "macroscopic" in the sense that, for 
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Figure 1. Asymptotic behavior of the mean number of real roots (N n ) of P n (x) in 
Eq. (J3j) as a function of a. These polynomials exhibit a condensation of their roots on 
the real axis for a > 2. 



large n, there is a finite fraction (N n )/n of the roots which are on the real axis. For 
L n {x) in Eq. ((2j) this fraction is exactly one. We thus say that there is a condensation 
of the roots on the real line, similar to a Bose-Einstein condensation where a finite 
fraction of the particles of a quantum-mechanical system (Bosons) condense into the 
lowest energy level. In the case of random polynomials, the roots play the role of the 
particles and the equivalent of the ground state is the real line. 

The purpose of this paper is to understand what types of polynomials lead to this 
condensation phenomenon. Of course, it is very difficult to address this question for 
any random coefficients a^. However, guided by the two examples above in Eq. (pQ) 
and Eq. (J2J, and in particular by the large k behavior of (o^), we introduce a family of 
random polynomials P n (x) indexed by a real a > defined by 

n 

P n (x) = J2 a kX k , (a 2 k ) = e~ ka , (3) 

k=0 

where are real independent Gaussian random variables of zero mean. While a = 
corresponds to Kac's polynomials, we recall that, for W n (x) in Eq. (pQ), (a\) oc e~ klogk 
and for L n {x) in Eq. ((2j), (a|) oc e - 2fc2l °s fc . Therefore, due to the extra logarithmic 
factor, these random polynomials are not exactly of the form introduced above ([3]). 
However, for a — > 1 + , one expects to recover the behavior of W n (x) in Eq. (pQ) while for 
a — > 2 + , one expects P n (x) to behave similarly to L n (x) in Eq. (j2J) : this is depicted 
schematically in Fig. [D 

Our main results can be summarized as follows. As a > is varied one finds 
three different phases. The first phase corresponds to < a < 1, where one finds 
that (N n ) ~ (2/7r)logn. In the second one, corresponding to 1 < a < 2, one has 



(N n ) ~ ^yj n a l 2 '. And in the third phase, for a > 2, one finds (N n ) ~ n. The 
condensation of the roots on the real axis thus happens for a > 2 and as one increases 
a, the condensation transition sets in at the critical value a c = 2. Furthermore, one 
finds that these real roots condense into a quasi-periodic structure such that there is, 
on average, one root in the interval [— x m+ i, — x m ) U [x m ,x m+ i], with x m = e% ma 1 , 
with 1 m < n. These different behaviors are summarized in Fig. [H By analogy 
with phase transitions of statistical systems the case < a < 1 can be considered as 
a high-temperature phase whereas a > 2 corresponds to the low-temperature (ordered) 
phase. Roughly speaking, one can consider our results as an interesting example where 
the transition from the high temperature where (N n ) oc logn (governed by a "a = 
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fixed point") to the low temperature phase where (N n ) oc n (governed by "a = oo" fixed 
point) happens through a marginal phase, for 1 < a < 2, where (N n ) ~ with an 
exponent <fi = a/2 which depends continuously on a. 

The paper is organized as follows. In section 2, we describe the general framework 
to compute the local density of real roots, which directly leads to (N n ). In section 3 
to 6 we then analyse separately the cases 0<a<l,a<2,a>2 and the "critical 
case" a = 2. In section 7, we give a qualitative argument to explain the condensation 
transition occurring at a c = 2 before we conclude in section 8. The Appendix contains 
some useful technical details. 



2. General framework 

First we notice that given that P n (x), as a function of x, is a Gaussian process, it is 
completely characterized by its two-point correlation function C n (x,y) 

n 

C n (x, y) = (P„(x)P„(y)> = ^ x " V k > ( 4 ) 

where we used the notation (...) to denote an average over the random variables a^. 
A central object involved in the calculation of (N n ) is p n {x), the mean density of real 
roots at point x. If we denote Aj, A2, X p the p real roots (if any) of P n (x), one has 
S(P n (x)) = Yli=i 3( x ~ Ai)/| J P^(Ai) I such that p n (x) can be written as 

Pn (x) = ][>(x - Ai)) = (\P^x)\6(P n (x)) 



1 

00 



dy\y\(5{P' n {x) - y)5{P n {x))) . (5) 

Under this form (jSJ), one observes that the computation of the mean density involves 
the joint distribution of the polynomial P n (x) and its derivative -P^(x) which is simply 
a bivariate Gaussian distribution. After Gaussian integration over y, one obtains 

_ y/Cn(x)(d n (x)/x + ggj - [c' n (x)] 2 
M) ~ 27TC n (x) ' ^ 



n 



-k a x 2k 



^n(^) C n (yX,X^j ^ ^ 6 
fc=0 

This formula ([6]) can be written in a very compact way [4] : 

Pn(^) = -\/ d u d v \ogC n {u,v) . (7) 

7T u=v=x 

Given that the random coefficients a& are drawn from a symmetric distribution, we can 
restrict our study of p n {x) on M + from which one obtains the mean number of real roots 

{N n ) as 

poo 

(N n ) = 2 / p n (x)dx . (8) 
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An important change of variable. We will see below that it is useful to consider 
these polynomials P n (x) in terms of another variable Y defined as 

/ 2 \ s=r 

Y = ( -logx J . (9) 

We denote p n {Y) the mean density of the real roots in terms of this new variable such 
that one has also (N n ) = J °° p n (Y)dY . For < a < 1 we will see that, for large n, 
most of the real roots of P n are located close to Y = n while for a > 1, the density 
extends over the whole interval Y G [l,n]. This change of variable ([HD is motivated by 
the following analysis. 

First we notice that C n (x, y) = Y^k=o e~ fcQ x fc ?/ fc in Eq. (pj is of the form C n (x, y) = 
c n (y/xy). Anticipating a saddle point analysis, one writes c n (x) as 

n n 

Cn (x) = e~ ka x 2k = ex P (-^(k, x )) , <P( k , x) = k a - 2k logx . (10) 

k=0 k=0 

Although 4>{k, x) is defined for integers k = 0, 1, 2, • • • , n, it is readily extended to the 
real axis and denoted <f)(u,x) = u a — 2m logx for u G IR + . The behavior of c n (x) is 
essentially governed by the behavior of <p(u,x) as a function of u (and fixed x). In 
particular, for a < 1, <j)(u,x) has a single maximum while for a > 1, it has a single 
minimum for u = u*(x) given by 

d n (j)(u*{x), x) = , d?<f){u*(x),x) = a(a - l)u*(x) a - 2 > , 



i 

ct-l 



u*(x) = ^— logxj . (11) 

The new variable Y introduced above in Eq. ([9]) is thus precisely Y = u*(x). As a 
consequence, the density behaves quite differently in both cases a < 1 and a > 1. 

For a < 1, most of the real roots on M + are located in [l,oo]. For fixed x > 1, 
4>(u, x) as a function of u in the interval [0, n] has a global minimum for u = n. Therefore, 
the sum entering in the expression of c n (x) in Eq. f fTDl will be dominated by the terms 
with k ~ n. The expansion of 4>(k, x) in Taylor series around k = n yields 

4>(k, x) = (f)(n, x) + {k — n) (an a ~ l — 2 log x) + • • • 

= (1 - a)n a - k(an a ~ l - 2 logx) + • • • , (12) 

where the higher order terms can be neglected in the large n limit because 
(9 J 0(n, x)/dui = 0{n a ~i) for j > 2. Thus, for a < 1 one has 



cjx) ~ e - {1 - a)n ' 



5>e-^-T , (13) 



k=0 



which, in terms of the rescaled variable x = xe~t nQ 1 j is the correlator of Kac's 
polynomials. From this observation ([TBI , one can straightforwardly obtain the mean 
number of real roots (N n ), this will be done in section 3. 

For a > 1, the situation is quite different and in that case, (f>(u,x) has a single 
minimum for u = u*(x) = (^logi)" 31 (Tni ). Besides, we will see below that the 
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main contribution to (N n ) on IR + comes from the interval 1 < x < exp (§^ a_1 ) where 
1 < u*(x) < n. In that case the sum entering in the definition of c n (x) in Eq. (fTUl) is 
indeed dominated by k ~ u*(x) and c n (x) can be evaluated by a saddle point calculation. 
For this purpose, one obtains after some algebra explained in the Appendix, a convenient 
expression of p n (x) as 



IX 



1 

7TX 



Z n k=0 (k-u*(x)ye 



-(f>(k,x) 



k=0 ' 



-<j>(k,x) 



J2 n k =o(k-u*(x))e 



-4>{k,x) 



-i 2 N 



fc=0 



-<fi(k,x) 



,(14) 



which is the starting point of our analysis for a > 1. For 1 < x < exp (|n a_1 ), one has 
u*(x) < n so that the sums over k in Eq. (fT4l ) are dominated by k ~ it*(x). The Taylor 
expansion of x) around this minimum reads 



(j)(k,x) = (f)(u*(x),x)+^2 

J'=2 



a(a — l)...(a — j + 1) 



(k-u*(x)) j lu*(x)] a -J .(15) 



For large x, u*(x) oc (logx) 1 ^ 0-1 - 1 is also large so that, to leading order in x, one can 
retain only the term corresponding to j '• = 2 in the Taylor expansion in Eq. (jTHI) . This 
yields, for large x 



g(k - u*(x)) exp (-<p(k, x)) 



(16) 



fc=0 



J^5((A; - u*(x))exp 



k=0 



■(k-u*(x))[u*(x)} 



a-2 



with g(z) = z or g(z) = 2 2 as in Eq. ( 1T41) . For later purpose it is useful to write 
u*(x) = [u*(x)\ + b with < b < 1, where L M *( a; )J is the largest integer smaller than 
u*(x) (i.e. the floor function). Performing the change of variable m = k — [u*(x)\ in 
the discrete sum ([TBI) , such that k — u*(x) = m — b one obtains the useful expression 

n 

^2g{k-u*(x))exp(-(f>(k,x)) (17) 



k=0 



n— [u* (x)] 

z -<t>(u*(x),x) g(m-b)exp 

m=— \u* {x)\ 



a(a — 1) 



{m-bf[u*{x)} 



a-2 



One clearly sees in expression f fT8l) that the behavior of this discrete sum, due to the 
term {u*(x)] a ~ 2 oc (loga;)^ -2 ^" -1 ) , will depend on the sign of a — 2. We will thus treat 
the three cases l<a<2, a > 2 and a = 2 separately. This will be done in section 4, 
5 and 6 respectively. 



3. The case < a < 1 



In that case, from the expression for c n (x) in Eq. ([TBI) , we can use the results of 
Kac's polynomials to obtain that most of the real roots will be such that, for large n, 
0{n~ l ) [TJ]. In other words, the real roots are distributed in a region 



_a a-l 

xe 2 — 1 
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of width 1/n around e? "° = 1 + ^n a ~ l + 0{n a ~ 2 ) and this distribution is exactly the 
same as the one for Kac's polynomials (corresponding to a = 0). The number of real 
roots is thus also the same and given by 

(AQ~-logn, (18) 

7T 

independently of a < 1. 
4. The case 1 < a < 2 

In that case [u*(x)] a ~ 2 — > for large u*(x) and one thus sees on the asymptotic 
expression in Eq. ([TBI that the discrete sum can be replaced by an integral. This 
yields, for large n and large x with x < exp (| n a_1 ) 

^g{k- u*{x)) exp (-0(Jfe, x)) ~ e"*^*^ / ff(y)e~ iii VV ,t * (x > a ~ a cfy .(19) 
fc=o 

Note that the prefactor e - ^"*^' 3 ^ is unimportant for the computation of p n (x) because 
it disappears between the numerator and the denominator in Eq. (fT4l ) and it will be 

omitted below. In particular, setting g(z) = 1 in Eq. (fT9l) one has 

i 

\2-Q!" 



exp (—<p(k, x)) oc v^r 

fc=0 



_a(a — 1) 

and similarly, setting g(z) = z 2 in Eq. ( TT9l l one has 



(20) 



a(a — 1) 



(21) 



yj(A: — u*(x)) 2 exp (—<fi(k, x)) oc 

fc=0 

while 5^fc =0 (^ — w*(x)) exp (— 0(fc, x)) ~ to lowest order in n. Therefore using the exact 
expression given in Eq. (fT4l) together with the asymptotic behaviors given in Eq. ( 1201 
|2T1) . one obtains the large x behavior of p n (x) as 

1 1 (2 \ 

Pn{x) -logx . (22) 

vrx ^/ a (a - 1) V" / 

For a clear comparison with the case a > 2 (which will be analysed in the next section), 

i 

it is convenient to write the density p n (Y), in terms of the variable Y = (^logx) a_1 , 
which reads, for 1 <C Y < n 



MY) ~ y^ZllY-i^ , (23) 

and in Fig. [2] a), we show a sketch of this asymptotic behavior ( 1231) of p n (Y) for 
1<7 <n. 

We can now compute (iV n ) = p n (x)dx. First, one notices that for a > 1, the 
series entering in the definition of c n (x) in Eq. (flOl) has an infinite radius of convergence 
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so that one readily obtains that f^p n (x)dx is of order 0(1) in the limit n 
Besides, for large x 3> et n< * \ one has (see also Ref. [7]) 

1 



X-i> 



a ~,a — 1 
C2 n 



TCX' 



OO. 



(24) 



(<>»> 71x2 

which implies that f °S n a-i p n {x) dx is also of order 0(1) in the limit n — ► oo. From these 
properties, it follows that the main contributions to (N n ) on IR + comes from the interval 
[1, e% n " X ] where the asymptotic behavior of p n (x) is given in Eq. ( |22i ). Therefore one 



has 



a „ot— 1 



7T 



a: 



a/2 



a: 



(25) 



where the factor 2 comes from the additional contribution coming from [— e^ n , — 1]. 
We thus have here an algebraic growth (N n ) oc n a l 2 with a continuously varying 
exponent a/2. This exponent tends to 1/2 as a — > 1 + , which is expected from the 
analysis of Weyl polynomials W n (x) in Eq. JTJ) for which {a\) oc e~ klogk (although the 

. Besides, from Eq. (i25l) . one also 



variance is not exactly of the form (al) = e~ k ' 
obtains that the amplitude of this term proportional to n a l 2 vanishes when a — > 1. We 
recall that for a < 1, one has instead (N n ) oc (^)logn ( fTST ). characteristic for Kac's 
polynomials. This suggests that this limit a — >• 1 is rather singular in the sense that the 
asymptotic behavior of (iV n ) for large n changes "discontinuously" from logn to y/n. 



5. The case a > 2 



In that case, the behavior of the discrete sum in Eq. (1T8|) . which enters in the 
computation of p n (x) ( fi~4l) is quite different. Indeed, in that case [u*(x)] a ~ 2 oc 
(logx)( Q ~ 2 ' ) / <Q ~ 1 ) — > oo for large x and therefore the leading term for large x in Eq. (TTHl) 
corresponds torn = if 6 < 1/2 or m = 1 in b > 1/2. Keeping these leading 
contributions, one has 



9 ( k - u*(x)) exp (-(j)(k, x)) oc g{-b) exp 
a(a — 1 



k=0 

+ ^(1 — b) exp 



a(a — r 



b 2 u*(x) a ~ 2 



■(l-b) 2 u*(x) a - 2 



(26) 



where, again, we have omitted the unimportant prefactor e ^ u *^' x \ Using this large x 
expansion (l26l) . one obtains p n (x) in Eq. (TT4l) as 



Pn{x) 



(irx) cosh 



a ("- 1 ) ya-2 



Y 



, ■ (1 - 26) 
In terms of the variable Y, the density p n (^) reads, 

a(a - l)r a " 2 



a 



logx 



i 

a-l 



(27) 



p n (F= L^j+6) 



27r cosh 



«(«-!) y a „2 



;i - 26) 



(28) 
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In Fig. [2] c), one shows a sketch of p n (X) for large Y < n given by Eq. (1251) : it is 
qualitatively very different from the case a < 2 (see Fig. [2ja)). Indeed, p n (Y) exhibits 
peaks centered around k + | for large integers 1 <C k < n. The height of these peaks is 
given by a(a — l)k a ~ 2 / (2n) whereas its width scales like k 2 ~ a . 

From p n (x), one can now compute the mean number of real roots. As in the case 
a < 2 (see Eq. (1231) and above), one can show that the main contribution to (N n ) comes 
from the intervals [-e^ n " , —1] and [l,e^ nQ ]. One thus has from Eq. (1251) 



(N n ) = 2 / p n (x) dx~2 p n {Y) dY (29) 
Jo Jo 

^ f 1 a(a-l)k^ n ^ dz 
~ > / f T do ~ > / — , 

f^Jo vrcosh ^Hk a - 2 (l-2b) ^J-M^i) fc «-2 7rcoshz 

and finally 

(JV n ) ~ n , (30) 

where we have used f_ dz/ coshz = it. This condensation of the roots on the real 
axis, characterized by the fact that (N n ) ~ n thus occurs via the formation of this 
quasi-periodic structure (see Fig. [2jc)). More precisely, this computation in Eq. ([291 
shows that for large fc, 2 f* +1 p n (Y) dY ~ 1 which means, going back to the original 
variable x, that there is, on average, one root in the interval [—Xk+i, —Xk] U [xk, Xk+i], 
with Xk = ea' 



a foot — 1 



6. The special case a = 2 



In view of the previous analysis, it is tempting to consider the fraction of real roots 
$ = lim n ^oo ( N n ) jn as an "order paramater". For a < 2, one has $ = whereas $ = 1 
for a > 2. One can however interpolate smoothly between these two limiting cases by 
considering the case a = 2 and introducing an additional real parameter p such that 

„2\ _ 



(31) 



Performing the same algebra as explained in the Appendix, one obtains the same formula 
as given in Eq. (1T31) with u*(x) = p^ 1 logx. The new variable is thus here Y = p" 1 \ogx 
and, setting Y = \ Y\ +b it is easy to see that the density p n (Y) is given by for 1 <^ Y < n 

1/2 



Pn(Y) 



P 
TV 



E 



rn = 



E 



oo 



-li(m—b) 2 



E 



rn= 



(m — b)e 



-/i(m—b) 2 



E 



oo 
rn = 



-li(m—b) 2 



(32) 



which is thus 1-periodic for all p. In Fig. [2jc), one shows a sketch of p(Y) for a = 2 
given by Eq. (l32l . For p — > 0, the density is almost constant and p n (^) ~ ^ x \J pj^ 



and the modulation of the density increases with p. For large /i, the sum in Eq. (1321) is 
dominated by the terms corresponding to m = and m = 1 and p„(y) is thus given by 
a formula similar to Eq. ( 1251 ) setting a = 2 and replacing y a ~ 2 by /i. For the average 
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Figure 2. a) : Sketch of p n (Y) (in arbitrary units) given in Eq. lj23| as a function 
of Y for 1 <C Y < n for a < 2. b) : Sketch of p n (Y) (in arbitrary units) given in Eq. 
(|32"1) as a function of Y for 1 <C Y < n for a = 2. c) : Sketch of p n (Y) (in arbitrary 
units) given in Eq. lj28j) as a function of F for 1 <C Y < n for a > 2. Here fc denotes 
an integer with 1 <g; k < n. 



number of real roots one has 



(JV n > oc \ it (33) 
n , jU>1 , 

which shows that this family of real random polynomials f l3~Tl) interpolate smoothly 
between the cases a < 2 ([25]) and a > 2 (l29l). 



7. A qualitative argument for the transition at a = 2 



This condensation of the roots on the real axis can be qualitatively understood if one 
considers the random polynomials (for x > 0) P n {X) = P n { x ) °f the variable Y, which 
one writes as 



--(k a - akY a ~ l ) 



(34) 



P n {Y) = Y,a k w(k,Y) , w{k,Y) = exp 

k=0 

and Sfe are i.i.d. Gaussian variables of unit variance. It is easy to see that the weights 
w(k, Y), as a function of k, have a single maximum for k = Y where the second derivative 
is proportional to k a ~ 2 . Thus for a > 2, the weights get more and more peaked around 
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this maximum for large k, whereas is typically of order 0(1). Therefore, given a 
large integer m, P n (fn) is, for a > 2, dominated by a single term corresponding to 
k = m. Consequently, the sign of P n {m) is essentially the sign of a m . This in turn 
implies that, if a m and a m +i have an opposite sign, P n (x) has, with a probability close 
to 1, a root in the interval [e^™" 1 , e^™" 1-1 ^ ]. In the case where a m and a m+ i have the 
same sign, the same argument shows that P n (x) has, with a probability close to 1, a root 
in the interval [— et( m+1 ) Q 1 ; _ e f( m ) a ]. One thus recovers qualitatively the result we 
had found from the computation of p n (Y) in Eq. (1291 ) where we have shown that P n (x) 
has, on average, one root in the interval [— e% ( m+1 ) Q 1 ; _ e f ( m ) a u [ e f ma 1 ; e f (m+i) Q x j_ 
This shows finally that P n {x) has, on average, (N n ) oc nreal roots. 

We also point out that our argument explains in a rather intuitive way the result 
obtained by Littlewood and Offord [H] for the random polynomials L n (x) ([2]). For 
these specific polynomials, defining xq = 0, x m = m m m\, they rigorously proved, using 
algebraic (and rather cumbersome) methods, that L n (x) has a root either on [x m ,x m+ i] 
if e m e m+ i = —1 or in [— x m+ x, —x m ] if e m e m+ i = 1. Our argument gives some insight on 
their intriguing result and allows to understand it in a rather simple way. 

8. Conclusion 

To conclude we have introduced a new family of random polynomials (J3j) , indexed by 
a real a. For these random polynomials, we have computed the mean density of real 
roots p n (x) from which we computed the mean number of real roots (N n ) for large 
n. We have shown that, while for < a < 1, (N n ) ~ (^)logn, the behavior of 
(iV n ) for a > 1 deviates significantly from the logarithmic behavior characteristic for 
Kac's polynomials. For 1 < a < 2, we have shown that (N n ) oc n a ^ 2 whereas for 
a > 2, (N n ) ~ n. This family of real random polynomials thus displays an interesting 
condensation phenomenon of their roots on the real axis, which is accompanied by an 
ordering of the roots in a quasi periodic structure : this is depicted in Fig. [2j 

Of course, the occurrence of this transition raises several interesting questions like 
the behavior of the variance of the number of real roots for large n as a is varied. It 
would be also interesting to compute the two-point correlation function of the real roots, 
which is a rather natural tool to characterize this periodic structure we have found. In 
view of this, we hope that this interesting phenomenon will stimulate further research 
on random polynomials. 

Appendix A. A useful expression for the mean density p n (x) 

In this appendix, we derive the expression for the mean density p n (x) as given in Eq. 
( 1T41) starting from Eq. ([7]). We first write c n (x) = (P n (x)P n (x)} as 

n 

Cn (x) = e -*(«*(*).«) e~ kk,x) , (A.l) 

k=0 
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where u*(x) is the location of the minimum of (j)(u,x) given in Eq. ([TT 



u*(x) = Qloga^ a 1 , (A.2) 

and 

(j)(u*(x),x) = (1 - a)u*(x) a (A.3) 

4>(k,x) =(f>(k,x)-(f)(u*(x),x) = fc a -a%*(x)] a-1 + (a-l)[u*(a:)] a . 

The correlator C n (x,y) = c n (y/xy) is given by Eq. (1A.1I) together with Eq. (1A.3I) where 
x is replaced by y/xy. All the dependence of C n (x, y) in x, y is thus contained in u*(y/xy) 
only. From its definition in Eq. (IA.2I) one has immediately 

d x u*(^5y) = 1 1 v -[u*(V^)] a " a , (A-4) 
a (a — 1J a; 

from which we obtain a set of useful relations 

%, y <f>(u*(Vw), VW) = -^TT) ^[^(v^)] 2 ^ (A.5) 

d x 4>(k, y/xy) = -{u*(yfxy) - k) 

d 2 x J(k,^xy) = 1 -\u\^)f- 2 . 
,y a{a — 1) xy 

For the computation of p n (x) from Eq. (JTJ) , it is useful to introduce the notation, for 
any function g{k) 

{g{k))z = E^Wexp(-feV5y)) (A 6) 

ELo ex P (-^(^^ V%y)) 

From C n (x,y) = c n (^fxy) and c n (x) given in Eq. f lA.lh one obtains 

d x d y \ogC n (y/xy) = -dl y (j)(u*(y/xy), y/xy) - (d x <p(k, y/xy)d y (j>(k, yfxy)) z 

- (dj{k, y/xy)) z (dj(k, y/w)) z - (dlj(k, Jxy)) z . (A.7) 

From the above relations in Eq. (|A.5|) , it is readily seen that the first and the last term 
in Eq. (I A. 71) cancel each other. Using the relation in Eq. j7j), one finally obtains the 
relation given in the text in Eq. (fl4l) . 
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